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Tap-proof long-distance quantum communication requires a deep understanding of the strong 
losses in transmission channels. Here we provide a rigorous treatment of the effects of beam wan- 
dering, one of the leading disturbances in atmospheric channels, on the quantum properties of light. 
From first principles we derive the probability distribution of the beam transmissivity, with the 
aim to completely characterize the quantum state of light. It turns out that beam wandering may 
preserve nonclassical effects, such as entanglement, quadrature and photon number squeezing, much 
better than a standard attenuating channel of the same losses. 
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Introduction. - The quantum properties of light open 
challenging perspectives for secure data transfer [Tj. 
Quantum communication is already used commercially, 
but its application is still restricted by relatively small 
distances. The main limitations are the strong losses over 
long distances. For example, in optical fibers typical val- 
ues are 0.2dB/km allowing us to preserve quantum 
entanglement over 100 km [3 . However, even this dis- 
tance is insufficient for global quantum communication. 
Some possibilities to improve this situation are based on 
protocols, which are rather stable against losses U, or on 
the idea of quantum repeaters [5]. 

A promising alternative is based on free-space channels 
for satellite-based global quantum communication, which 
can bridge any distance on the Earth. In recent exper- 
iments quantum entanglement of light has been demon- 
strated, after transmission through a free-space channel 
of 144 km with atmospheric losses of 32 dB [W, T. Atmo- 
spheric losses are instable and permanently fluctuating. 
The related noise is expected to substantially diminish 
the quantum properties of the transmitted light. How- 
ever, such fluctuations may even be useful to preserve 
quantum entanglement [8] . 

Noise effects in free-space channels have been exten- 
sively studied [HI HH] i including models with random fluc- 
tuations and stochastic modulation of the intensity |111 - 
[14] . More recently, the method of fluctuating-loss chan- 
nels has been introduced (TS] , based on fluctuations of the 
complex transmission coefficient, whose absolute value is 
bounded by one [121 HZl • Beyond the earlier models [llr 
[13], the latter work in Ref. |15] yields a consistent de- 
scription of quantum channels. 

The probability distribution of the transmission coef- 
ficient (PDTC) is the main characteristic of fluctuating- 
loss channels. However, the lack of a consistent model for 
the PDTC prevents us from exploiting the advantages of 
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this theory under realistic conditions. The phenomeno- 
logical log-normal model was studied for the random in- 
tensity modulation jl3) . but it incorrectly describes the 
distribution tail with the transmission coefficient being 
close to 1. This domain, however, is most important for 
preserving the quantum effects of light needed for tap- 
proof communication. 

For a realistic description of the quantum effects of 
light in free-space channels, the derivation of a PDTC 
from first principles is indispensable. Among the un- 
wanted disturbances in free-space channels we will con- 
centrate on the important phenomenon of beam wan- 
dering, caused by turbulence [HI [TO] and unstable adjust- 
ment of the radiation source. We obtain an analytical ex- 
pression for the transmission coefficient of an aperture. 
Based on a normal distribution of the beam position- 
ing |14| , we derive an analytical form for the correspond- 
ing PDTC. Later on, we study the effects of beam wan- 
dering on quantum phenomena, such as photon-number 
and quadrature squeezing, and violations of Bell inequali- 
ties. We show that despite of the extra noise, fluctuating 
loss channels may preserve such nonclassical properties 
much better than similar constant-loss channels. This 
result gives a strong argument to favor free-space quan- 
tum communication. 

Fluctuating-loss channels.- The description of losses 
in linear quantum optics usually connects annihilation 
operators of the input and output fields, ain and Oout, 
by the standard input-output relation, Oout = T^Oin + 
^yi — T'^c cf. [TO], Appendix A. The operator c describes 
modes of the environment being in the vacuum state. 
For our further considerations we assume the absence 
of dephasing. This applies to many experiments, such 
as photodetection [TO] and polarization analysis for test- 
ing Bell inequalities f7^. One may also perform homo- 
dyning, using a local oscillator copropagating with the 
signal [TO] [10]. Hence, the transmission coefficient T is 
real and positive. Additionally, preserving the commu- 
tation relations implies that Te [0, 1]. In fluctuating-loss 
channels, T is a random variable. The corresponding 
quantum-state input-output relation attains its simplest 
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form in the Glauber-Sudarshan P representation 
(cf. [18], Appendix B), 
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where Pin (a) and Pout (a) are the input and output 
P functions, respectively. Fluctuations of losses are de- 
scribed by the PDTC, V{T), which we will derive for the 
scenario of beam wandering. 

Aperture transmission coefficient.- For weak absorp- 
tion, beam- wandering losses are dominant. Typically 
they are caused by aperture truncation of the light at 
the receiver. This implies that T in Eq. ([T]) is the aper- 
ture transmission coefhcient of a light pulse. We consider 
the pulse as a superposition of Gaussian beams |21j with 
different wave numbers k, propagating along the Z axis 
onto the aperture plane at distance Zap from the source. 
If beam deflection is caused by imperfect adjustment of 
the radiation source or by weak turbulence, the beam in- 
cidence is normal to the aperture plane to a good approx- 
imation. The beam center is deflected by the distance r 
from the aperture center, see the inset in Fig. [l] 




FIG. 1. (color online) Transmission efficiency, T^, vs the 
beam-deflection distance, r, for difl^erent values of the beam- 
spot radius W. Solid and dashed lines correspond to the 
numerically calculated integral ([3| and the analytical approx- 
imation Q, respectively. The inset schematically shows the 
aperture of radius a and the beam-spot of radius W. 

The transmission efficiency (square of the transmission 
coefficient) of a Gaussian beam reads as 



dx dy\U{x,y,z^p;k)\" 
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where A is the area of aperture opening, U{x,y, Zi^p] k) 
is the Gaussian beam normalized in the XY plane. In 
practice the transmission coefficient T{k) approximately 
equals the pulse transmission coefficient in Eq. ([T]). That 
is, the pulse spectrum is irrelevant, T « T(fco), with fco 
being the carrier wave-number. For a Gaussian beam 
with the spot radius W at the aperture plane the trans- 
mission efficiency is given by the incomplete Weber inte- 
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a is the aperture radius and I„ the modified Bessel func- 
tion. Details can be found in 18j, Appendix C. 

The integral ^ can be calculated numerically [14j . 
However, for the evaluation of the PDTC we propose an 
approximate analytical expression of the form (cf. |18j . 
Appendix D) 



= exp 
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where Tq is the maximal transmission coefficient for the 
given beam-spot radius W] A and R are the shape and 
scale parameter, respectively. The constants Tq, A, and R 
are obtained from the following procedure. We consider 
the transmission efficiency as a function T'^(r) of the 
beam-deflection distance r. For the integral (Is]) one can 
calculate analytically the particular values T^fi), T'^{a), 
and {dT^{r)/dr)r=a cf. [22] ■ From the condition T^{0) = 



Tg one gets 



Tq^ = 1 - exp 
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The shape and scale parameters, A and R, are obtained 
from two algebraic equations for the values of and its 
derivative at the point a. 
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In Fig. [T] we compare the numerical result for T'^[r) with 
the approximation (|4])-([7]). The maximal relative mean 
quadratic error appears for W=Q.2ia and is 1.85%. This 
gives a reasonable accuracy for the applied technique. 

Probability distribution of the transmission 
coefficient.- We suppose that the beam-center po- 
sition is normally distributed with variance cr^ around a 
point at the distance d from the aperture center. In the 
case of beam wandering caused by imperfect adjustment 
of the radiation source this variance is cr^ "'S-^^^j 
where is the variance of the source-deflection angle. 
Likewise, in the case of weak atmospheric turbulence, 
cr2 w 1.919 C2z3(2Wo)-i/3 [10], where Cl is the index- 
of-refraction structure constant, Wq is the beam-spot 
radius at the radiation source; for more recent results 
see [23]. Generally, both variances should be added up. 

Based on the above assumption, the beam-deflection 
distance r fluctuates according to the Rice distribu- 
tion [TT] with the parameters d and a. Provided that the 
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transmission coefficient T is approximated by Eq. Q, 
the PDTC is given by the log-negative generalized Rice 
distribution, 
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for T G [0,Tq], and 7'(r) = else. In the particular case 
when the beam fluctuates around the aperture center, 
d=0, this distribution reduces to the log-negative WeibuU 
distribution. 
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for T e [0,To], and V{T) = else. More details about 
properties of the PDTC can be found in [TS], Appendix E. 

The PDTC ^ and (|9| may have singularities at the 
points r=0 and T—Tq. Evidently, in any experimental 
reconstructions of the PDTC this singularity cannot be 
observed directly due to finite values of the sampling- 
data number and the sampling interval for T. For this 
reason, we will use the exceedance (tail distribution). 
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which is the probability that the transmission coefficient 
exceeds the value T. Plots of this function are presented 
m Fig. [2) In the following we will show that contributions 
from tails with large T are important for preserving non- 
classical properties of the transmitted light. This implies 
that the exceedance may also characterize the feasibility 
of quantum protocols. 
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FIG. 2. (color online) The exceedance (tail distribution) 

cf. Eq. (10 1, is shown in the logarithmic scale for different 



values of the beam-spot radius W, and standard deviation of 
beam deflection a. Solid and dashed lines correspond to d=0 
and rf=cr, respectively. 



Nonclassical effects.- Input-output relation ([T]) and 
PDTC ([s]), ([9]) enable to consistently characterize quan- 
tum states of the transmitted light. We will use this pos- 
sibility for studying disturbances of nonclassical effects 
by beam wandering. In this context the fluctuating-loss 
channel with the mean transmission efhciency (r^) is 
compared with a standard attenuating channel with the 
transmission efficiency T'^= (T'^y 

A remarkable example is the experimental violation 
of Bell inequalities, after transmission over a 144 km at- 
mospheric channel with a mean loss of 32 dB |7] • As it 
has been reported in Ref. [5], this result is very sensi- 
tive to stray-light and dark counts. For example, such 
a violation is impossible for constant attenuation of the 
same size and the mean number of noise counts of 10~^ 
(instead of 0.5x10^^ in Ref. 0) cf. the dashed line in 
Fig. |3] Indeed, in this case the required simultaneous 
clicks of the detectors are mainly caused by stray-light 
and dark counts. For fluctuating losses events with large 
T, i.e. in the tail of the PDTC, ensure the simultaneous 
clicks of detectors from an entangled-photon source. A 
problem is that the earlier models, e.g. the log-normal 
PDTC, do not give a correct behavior for the tail. 
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FIG. 3. (color online) The Bell parameter is shown in de- 
pendence of the squeezing parameter of a parametric down- 
conversion source. The sohd line shows the situation for beam 
wandering, with a standard deviation a = 28.5a of beam de- 
flection and a beam-spot radius of = 1.1a, leading to a 
mean loss of 32 dB. The receiver-module losses are supposed 
to be 9dB. The dashed line shows the behavior for constant 
losses of the same size. The mean number of stray- light and 
dark counts is assumed to be 10~^. The shaded area indicates 
the violation of the Bell inequality. 

Hence it is important to verify that the tail of a con- 
sistent PDTC really improves violations of Bell inequal- 
ities, compared with the case of constant losses. In the 
present consideration we suppose that (on average) 32 dB 
of fluctuating losses are caused by beam wandering only. 
In fact, other phenomena significantly contribute to the 
measured statistics as well 6j. Here we effectively model 
all of them by the PDTC We also include into con- 
sideration 9 dB losses of the receiver module. 

We provide calculations similar to Ref. [8 but for the 
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PDTC ^ cf. also IS], Appendix F. The result is given in 
Fig. [3] It clearly shows that the entanglement can survive 
in case of beam wandering. Under similar conditions it 
disappears for constant losses of the same size. This find- 
ing is of vital importance for long-distance quantum com- 
munication based on entanglement and discrete-variable 
coding. 

The need to observe clicks at two polarization ana- 
lyzers post-selects also events with large T in the Bell- 
inequality test. In continuous-variable quantum commu- 
nication based on squeezed states |24| , such a natural pro- 
cedure does not occur. In order to preserve quadrature 
squeezing we should additionally monitor the transmis- 
sion coefficient and post-select the events with T exceed- 
ing a certain value T,„in- Likewise, similar procedure can 
be applied for improving the photon-number squeezing. 

Monitoring the transmission coefficient can be per- 
formed with the technique proposed in Ref. |20j , where a 
test pulse (also used as a local oscillator [TO]) is sent in 
the mode orthogonally-polarized to the signal. By appro- 
priately choosing Tmin we can, in principle, preserve any 
value of squeezing. However, this is limited by the occur- 
rence of the corresponding events. Hence the feasibility 
of this procedure must be analyzed. 

To be definite, assume that the input light is 
amplitude-squeezed by 6 dB and coherently displaced by 
an amplitude of 10. Similar to Bell inequality test, we 
suppose that beam wandering leads to 32 dB mean losses. 
In standard attenuation channels such losses result in 
negligibly small squeezing: 8.8xlO^*dB for the quadra- 
ture and 8.6xlO~^dB for the photon number. In Fig. |4] 
we present the value of squeezing in dependence on the 
exceedance (10) for the point Tmin, for details of calcula- 
tions cf. .^jjpendix G. It is clearly seen from the plots 
that approximately one event from lO'' is characterized 
by the high level of squeezing. Hence, under such con- 
ditions, the technique is still feasible. This conclusion is 
of importance for long-distance quantum communication 
based on continuous-variable coding. 

Summary and Conclusions.- Starting from the 
input-output relations for general quantum states we con- 
sistently describe the quantum properties of light under- 
going beam wandering in free-space channels. From first 
principles we have analytically derived the main charac- 
teristics of such channels - the probability distribution 



of the aperture transmission coefficient. It is given by 
the log-negative generalized Rice distribution. When the 
centers of beam-wandering and aperture coincide, this 
function reduces to the log-negative WeibuU distribution. 
The proposed distribution does not include phenomena 
other than beam wandering, which also may occur in tur- 
bulent media. However, our approach is the first model 
of free-space channels which consistently applies to the 
quantum domain. We believe that further generalizations 
will be possible. 

It has been shown that fluctuating losses associated 
with beam wandering can be more favorable for trans- 
ferring quantum properties of light, compared with con- 
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FIG. 4. (color online) The reachable values of quadra- 
ture and photon-number squeezing are shown as a function of 
the corresponding exceedance. The input light is amplitude- 
squeezed by 6dB and coherently displaced by an amplitude 
of 10. The mean losses are 32 dB with the same parameters 
of beam wandering as in Fig. [Sj The dependence of Tmin on 
the exceedance is also shown. 



stant losses of the same size, e.g. in fibers. This includes 
phenomena such as entanglement, photon- number and 
quadrature squeezing. Their conservation plays a cru- 
cial role for free-space quantum communication based on 
discrete or continuous variables. 
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SUPPLEMENTAL MATERIAL 

Appendix A: Operator input-output relation 

The description of linear losses in the absence of de- 
phasing is usually based on the operator input-output 
relation, 



Ta„ 



(Al) 



where Sin and Oout are the input and output field- 
annihilation operators, respectively, c is the operator of 
environment modes, T e [0, 1] is the transmission coeffi- 
cient. Moreover, T is assumed to be a random variable 



in the case of fluctuating- loss channels. Expression (Al) 



can also be considered as the input-output relation in the 
Heisenberg picture of motion. Here we will show how to 
derive this relation from the first principles for free-space 



channels and to get an explicit form for the transmission 
coefficient. 

We suppose that the atmosphere absorption is negligi- 
ble in the channel. In this case the light is scattered, so 
that parts of it do not reach the receiver aperture, which 
leads to losses. The positive-frequency part of the field 
operators before the aperture (the input field), yli„(r, t), 
and after the aperture (the output field), ^out(r,i), sat- 
isfies the wave equation, 

(r) ^iin/out (r, t) - Aii„/out (r, t) = 0, (A2) 

where c is the speed of light in the vacuum and 

n(r) = 1 + 5n{v) (A3) 

is the atmosphere index of refraction. These fields are 
related to each other as 



+00 



^out(r,i) 



dVr(r',r) Ai„(r',t) + C(r,t), (A4) 



where T(r',r) is the aperture transmission function, 
C{y, t) is the noise operator associated with the field scat- 
tered by the aperture. 

Let the input field be initially prepared in the form 
of a pulse, which corresponds to the normalized non- 
monochromatic mode ?7in(r,t). After passing the aper- 
ture this pulse transforms to the form 



+ 00 

j dVT(r',r) [/i„(r',i). 



(A5) 



The normalization of this pulse gives the corresponding 
output nonnionochromatic mode. 



+ 00 

y"dVT(r',r) [/i„(r',i) 



C/out(r,0 = 

/ -\-oo +00 

j A\ j A\'T{t',v) f/i„(r',t) 

y— oo — oo 

. ' 

which is detected and analyzed at the receiver. The func- 
tions f/in/out(r, also satisfy the wave equation (A2). 



The operators Ai,i/out(r, t) can be expanded into series 
of the nonmonochromatic modes, 



^in/out(r,i) — ain/out^^in/out(r, i) 

+c=o 



(A7) 
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in) 



in/out in/out 



where U^^J^^^{r,t) form a basis of the orthogonal com- 
plement of a mode C/in/out(r, i): and a[^j^^^ are corre- 
sponding field- annihilation operators. Let us substitute 
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Eq. ( A7 1 in Eq. ( A4 ) . Utilizing the orthogonality condi- 
tion, 



+00 



(A8) 



mode normalization, and Eq. ( A6 ) , we get operator 



input-output relation (Al) with the transmission coef- 
ficient, 



T ■ 



I +00 +00 

1 d^r 1 d3r'T(r',r) [/i„(r',t) 



(A9) 



being the normalization constant for the output pulse. 



cf. Eqs. (A5) and (A6). All terms containing the field 
annihilation operators are included in the noise oper- 



is the corresponding characteristic function. By using 
Eq. (Al), this allows one to derive the relation between 



the characteristic function of the attenuated and the in- 
put field, $t(/3) and <I>in(/3), respectively, as 



*t(/3) = $in(T/3). 



(B3) 



Here we utilize the fact that the characteristic function 
of the environment modes is equal to unity since all these 
modes are assumed to be in the vacuum state. Substi- 



tution of Eq. (B3 



relation between t 



in Eq. (Bl) gives the input-output 



le corresponding P functions. 



1 



PT(«) = ^Pi„(|). (B4) 
Now we have to take into account the fact that the 
transmission coefficient T fluctuates randomly. Conse- 



ator c of the input-output relation (TaTI) . The correspond- quently, we have to average relations ([BSf and (|B4|) with 



ing modes also contribute to the output pulse, however, 
all of them are in the vacuum state. 



the probability distribution of the transmission coeffi- 
cient (PDTC), 7'(r). The resufiing input-output rela- 
tions read as 



Appendix B: Quantum-state input-output relations 

In this Appendix we discuss the transformation of the 
operator input-output relation (Al) to the Schrodinger 
picture of motion, to obtain the corresponding density 
operators. The easiest way to solve this problem is based 
on the Glauber-Sudarshan P representation. In this rep- 
resentation the quantum-state input-output relation re- 
sembles the corresponding relation in classical optics. 

Let us remind that the P function, is related to 

the corresponding density operator q as 

+00 

P(a) = -1 y d2/3$(/3)e"'5*-"*^ (Bl) 



where 



= Tr {q exp exp (-a/3*)] (B2) 



*out(/3) = / dr7'(T) $i„(T/3) 



(B5) 



dTP(r) —Pi, 



'(r) ' 



(B6) 



where Pout (a) and <I>out(/?) are the P function and the 
corresponding characteristic function, respectively, of the 
output field. 



The quantum-state input-output relation ( B6 ) can be 
easily rewritten in any other representations. For ex- 
ample, the corresponding Wigner functions. Win (a) and 
Wout(a), are related to each other as 



Wout(a) 



1 -1-00 
dTP(r) 



7rT2 (1 - T2) 



exp 



1 -T2 



(B7) 



which requires an additional /3 integration. It is also use- 
ful to have the corresponding relation for the normal- 
ordered moments, Af„ 



Tr 



Igdt^a") 



(B8) 



which directly leads to Eqs. (G6) and (G7) of Ap- 
pendix [G| 



Appendix C: Pulsed Gaussian beam 



In this Appendix we specify the form of the transmis- 
sion coefficient (A9) for the case of the input pulse pre- 



pared as a superposition of Gaussian beams, see Ref. 
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with different wave numbers fc, 



(CI) 



1 



+ CXD 



dkC{k) Sin{x, y, z; k) exp[ifc (z — ct)] 



where C{k) is the pulse spectrum, Sin{x,y, z;k) is the 
input-beam envelope. The beams propagation direction 
corresponds to the Z-axis. Many laser-based sources 
emit light, which can be approximately described by 
Eq. (CI). The non-normalized pulse after the aperture, 
cf. Eq. (A5), has the form similar to Eq. (CI I, 



+00 

J dVT(r',r) C/in(r',t) 



(C2) 



+ 00 



/27r 



dkC{k) S'out (2;, Vi z; k) exp[zfc (z — ct)] , 



where S'out(a;, y, z; fc) is the ou tput -beam envelope. 
Since the pulses (CI) and (C2) obey the wave equa- 



tion (A2), the input- and output-beam envelopes in the 



paraxial approximation satisfy the equation 

Sin/out{x,y,z;k) ^Q, (C3) 



for details see Ref. |2T]. Here A_l is the transverse 
part of the Laplace operator. This equation pre- 
serves the norm in the (X, Y) plane, i.e. the relation 



+00 

J dxdy\Si, 



i/out 



{x,y, z;k)\ does not depend on z and 



equals to unity for the input pulse. The statement 
that the source irradiates Gaussian beams means that 



Sin{x,y,z;k) is the solution of Eq. (C3), with the follow- 
ing condition at z 



0: 



Sin{x,y,0;k) 



(C4) 



7rW§ 



exp 



2W2 



2i?n 



where Wq, Rq, and ^'o are values at the transmitter for 
the beam-spot radius, curvature radius of the wavefront, 
and the Gouy phase, respectively. 



We substitute Eq. (C2) in Eq. (A9) and take into 
account the fact that Sout{x,y, z; k) is slowly varying 
with respect to z, comparing with exp[ifcz]. Hence, 
in the corresponding integration we can approximate 
Sont{x,y,z;k) by Sont[x,y,z^p\k), where Zap is the po- 
sition of the aperture plane at the Z-axis. This results 
in the fact that the transmission efhciency (square of the 
transmission coefficient) has the form of decomposition, 



+00 



dk\C{k)\^T^ik), 



(C5) 



of the beam transmission efficiencies, 

+00 



T^{k)= / da;d7;|S'out(a;, y, Zap;A:)| 



(C6) 



with different wave numbers fc. 

At the aperture plane, z = Zap, the input and output 
beam envelops are related to each other as 

c „ , _ ( Sin{x,y,z^p;k) , foT{x,y) e A 
^outia;,y,Zap,Kj - I for (a;, y) ^ ' 

(C7) 

where A is the area of aperture opening. This implies 
that Eq. ( C6 ) reduces to 



T'^{k) = I dxdy |5i„(x, y, Zap;fc)|^ 

J A 



(C8) 



Here the integration is taken over the area of aperture 
opening. 

In practice, the shape of the pulse spectrum, |C(fc)|^, is 
much narrower than the shape of the beam transmission 
efficiency, T^{k). Besides, as it follows from the normal- 
ization condition for the input pulse, 



-t-00 

j dfc|C(fc)|^ = 1. 



(C9) 



This yields that the pulse transmission efficiency (C5) 



is approximately equal to the beam transmission effi- 
ciency (C8), T^«T^(fco), where fco is the carrier wave- 



number. Remaining that the Gaussian-beam amplitude 
is 

U {x, y, z; fc) = Sin{x, y, z; fc) exp(ifcz) , (CIO) 
the transmission efficiency can be equivalently given by 



« r2(fco) 



dx d?/|C/(x, y, Zap;fco)| 



(Cll) 



cf. Eq. (2). 

Therefore, in order to find the transmission coefficient 
one has to resolve Eq. (C3) with the boundary condi- 
tion (C4). For Sn{r) =0 the solution is obtained by re- 



placing the parameters Wq, Rq, and ^'o in the boundary 



condition (C4) with the corresponding values W, R, and 



vj*, which depend on z, see Ref. [21j for details. We model 
the effect of i5n(r) in Eq. (C3) by deflection of the beam 
from the aperture center, which works for the week tur- 
bulence only. Moreover, the beam can also be deflected 
due to the imperfect source adjustment. Since the aper- 
ture form is circular, all directions appear to be equiva- 
lent. Hence, we can consider, without loss of generality, 
the beam-deflection by the distance r in the direction of 
X-axis. The obtained solution, 



Sin{x, y, z; fc) 

1^ 



(C12) 



exp 



2W^ 



ik 
2R 



{x - r) +y^ 
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should be substituted in Eqs. (ClOl and (Clll. This 



results in the transmission coefBcient expressed in terms 
of the incomplete Weber integral, which is approximately 
evaluated in Appendix [Pj 



rj-iZ 



-e ^ »^ 



dppe'^^lo['i^), (C13) 



approximate expressions for T'^ir) and their derivatives 
dr^(r)/dr at the point r—a. The corresponding integrals 
for the exact relation (D2) are evaluated explicitly as, 
cf. Ref. [221, 



T 



(D5) 



where a is the aperture radius, and Iq is the modified 
Bessel function, cf. Eq. (3). 



Appendix D: Incomplete Weber integral 

In this Appendix we give some technical details related 
to the approximate analytical evaluation of the integral 
in Eq. (CIS), see also Eq. (3). First note that the incom- 
plete Weber integral is defined as, see Ref. [22] . 



Qn{x, z) = (2x) 



-n— 1 



dtt 



n+l 



exp 



(Dl) 

where I„ is the modified Bessel function. In terms of this 



special function Eq. (CIS) can also be given by 

,.2 



rpZ 



^ r ra 
2 4 



(D2) 



An alternative representation for Qnix, z) consists in the 
Lommel series of the modified Bessel functions. However, 
this series converges slowly and for practical purposes one 
usually uses the numerical integration. 



We derive an analytical approximation for Eqs. (CIS) 



and (D2 1. For this purpose we consider the transmission 
efficiency, T^, as a function T^(r) of the variable r. This 
function can be effectively approximated as 



T2 (r)=r2exp 



R 



(DS) 



cf. Eq. (4), where the parameters Tq, i?, and A (maximal 
transmission coefficient, scale, and shape parameters, re- 
spectively) can be obtained with using some properties 
of the incomplete Weber integral. 

The maximal transmission coefficient Tq is obtained by 
com parin g the exact, cf. Eq. (CIS), and approximate, cf. 



Eq. (DS), expressions for T^(r) at the point r=0. In this 



case the integral ( C1S[ ) is evaluated analytically. Using 
the straightforward calculations, the maximal transmis- 
sion efficiency. 



T2 = 1 



exp 



(D4) 



cf. Eq. (5), is obtained explicitly. 

Similar consideration can be applied for scale, i?, and 
shape. A, parameters. In this case we compare exact and 



\ dr 



4a 



-4 



2a^ Aa^ 



V 



(D6) 



These expressions can be equated with the same quanti- 
ties evaluated from the approximate form ( DS ) , 



T2(a) = To^exp 



(D7) 



/ dr^(r) 

dr 



A /a\^-i 



R \R 



exp 



f- 
\R 



(D8) 

The obtained system of two algebraic equations for R 
and A yields 



A=J 



exp[-4,fi]li(4^) 
W^M-cxp[-4|^]Io(4^) 

Inf 



2Ti 



(D9) 



R^a 



In 



l-cxp[-4^]Io(4^) 



l-exp[-4^]Io(4f^, 



(DIO) 



cf. Eqs. (7) and (8). 



Appendix E: Beam-wandering PDTC 

In this Appendix we present some technical details re- 
lated to deriving the beam-wandering PDTC, cf. Eqs. (8) 
and (9), and discuss mathematical techniques for oper- 
ating with this object. The basic assumption is that the 
beam-center position, {xo,yo), randomly fluctuates ac- 
cording to the two-dimensional Gaussian distribution. 



27rcr2 



exp 



2a2 



(El) 

where a is the standard deviation of beam deflection and 
A={dx,dy) is the mean position of beam center. This 
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implies that the beam-deflection distance r=yjx'^ + 
fluctuates according to the Rice distribution, 



p(r;rf, a) 



rd 



exp 



2a2 



(E2) 



where d=^ + is the distance between aperture and 

fluctuation centers, Iq is the modified Bessel function. 
The random variables r and T are related to each other 



via Eq. (D3|, cf. also Eq. (4). This yields that the proba- 



bility distribution of the transmission coefficient (PDTC) 
reads as 



V{T)-- 



2^ ^91 ^" 



(Rd 



2 In 



To 
T 



X exp 



(E3) 



for T e [0,To], and V{T) = else, cf. Eq. (8). For 
this distribution we refer to as the log-negative gener- 
alized Rice distribution. This is also important to note 
that log- losses 9^ — In T'^ are distributed according to the 
generalized Rice distribution. 



v{e) 



X exp 



\-i iRd. 



l-^(R^{e-e,y+d^) 



(E4) 



for 6* > 6*0 and V{9)^Q else, where 6*0=- InTp^. A partic- 
ular case of Eqs. ^ and (E4| for d=0 is the log-negative 
Weibull distribution for the transmission coefficient, 



V{Ty 



exp 



2a 



(E5) 



for T e [0,ro], and V{T) = else, cf. Eq. (9), 
and the Weibull distribution with the shape parameter 

(-\/2(T/i?) and the scale parameter 2/A for log-losses. 



no) 



exp 



1 

'2^ 



R^O-0o) 



(E6) 



for e>eo and V{e)^0 else. 

Integrations over the PDTC ^ and ^ is convenient 
to perform with respect to the measure dr{T). Here 



r(T) = i? 21n 



To 
T 



(E7) 



inverts Eq. (D3). In practice this means that the mean 



value for the function /(T) of the transmission coefficient 
T can be obtained as 

1 +00 

dTV{T)f[T]^ J drp{r;d,a)f[T{r)], (E8) 



where p{r; d, a) is the Rice distribution. 



Fluctuations of the transmission coefficient can also be 
characterized by the cumulative probability distribution. 



1 

T{T) = j AT'V{T') . 



(E9) 



For our purposes, however, it is convenient to use the 
exceedance (tail distribution), which is defined as 



± 

T{T) = 1 - j'(r) = J dT' V{T') , 



(ElO) 



cf. Eq. (10). Clearly, the exceedance is equal to the 
probability that the transmission coefficient exceeds the 
value of T. 

The exceedance of log-negative Rice distribution (|8| i s 
expressed in terms of the incomplete Weber integral (|D1|) , 



T{T) = 1- exp 



Qc 



2ct2' 



iT)d\ 



(Ell) 



where r{T) is given by Eq. (E7). For this special function 
we can apply the technique presented in Appendix [D] 
This results in an approximate expression. 



TiT)^l-To{T) exp 



where 



MT) = 1-exp 



D{T)J 



tiT) 



2cr2 



(E12) 



(E13) 



D{T) = r(T) { In 



2MT) 



l-Io(^)exp[-^]J 



(E14) 



'■(T) 



l-Io(^)exph^] 



X < In 



2-Fo(T) 



l-Io(^)exp[-^] 



(E15) 



For d—^ the approximate expression ( E12[ ) reduces to 

r\T) 



^(T) = l- J-o(r)==exp 



2a2 



(E16) 



which also coincides with the exact form of the ex- 
ceedance for the log- negative Weibull distribution ([9]). 
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Appendix F: Bell-inequality test 

In this Appendix we remind technical details from 
Ref. [8] concerning calculations of the Bell parameter 
for the light generated by a parametric down-conversion 
source, transmitted through the atmosphere, and then 
treated by polarization analyzers with on/ofF detectors. 
We consider experimental configuration used in Ref. [7], 
cf. also Fig. |F1| In this scenario both receivers are situ- 
ated at the same place and photons are split by a small 
time interval. 



Receiver 

Transmitter Telescope pgg 

HWP 





Delay Line 

FIG. Fl. Setup for checking Bell inequalities with the 
light transmitted through a 144 km atmospheric channel, 
cf. Ref. [7]. A parametric down-conversion source produces 
entangled photon pairs. These photons, separated by a small 
time interval, are sent to the receiver. After collecting by the 
telescope, photons are split by a beam-splitter BS and then 
treated by the corresponding polarization analyzers. Each 
polarization analyzer consists of (i) half-wave plate, HWP, 
which changes the polarization angles 9a and 6b; (ii) polar- 
izing beam-splitter, PBS; (iii) on/off detectors Dt^, Dtq for 
the transmitted light, and Dr^, Dr^ for the reflected light. 

The Bell inequality in the Clauser-Horne-Shimony- 
Holt (CHSH) form states that the parameter 



-E 



i(l) /,(2) 



(Fl) 



also referred to as the Bell parameter, cannot exceed the 
value of 2 in any local realistic theory. Quantum theory 
violates the Bell inequality. The correlation coefficients 
E{0a,O-b) in Eq. (IfI]) are defined as 



where 



,03) -Pd iffercnt {0a,9b) 



fforont {OA, Ob 



(F2) 



Psamc(eA,0B) = ^Ta ,Tb (^A , ) + Pr^.Rb (^A, ^b) (F3) 

is the probability to get clicks on detectors in the same 
channels of the polarization analyzers, and 



different \ 



i)=P 



TA,RBlfA,CB 



)+Pl 



Ra,Tb 



9a,0b) 
(F4) 

is the probability to get clicks on detectors in differ- 
ent channels. Here ^1^,153 (6'a, ^b) is the probability of 
registering clicks at the detectors za = {Ta,Ra} and 



iB = {Tb,Rb} for the polarization angles 9a and 6*3. 
According to the photodetection theory, these probabili- 
ties are given by 



-t-00 

0b) = E Tr(ll(^)ll(:)ll(>fJ^), (F5) 

n,m— 1 



^A 7^ JA, *B 7^ JB, where g is the density operator, 

n(") ('^"^AtB) +^)' 



«A(B) 



exp [-?7n»A(B) - ^] ■ (F6) 



is the positive operator-valued measure for the detector 
*A(B): N is the mean number of stray-light and dark 
counts [Semenov et al., Phys. Rev. A 78, 055803 (2008)], 
77 is the detection efficiency, : • • • : means normal ordering, 
and 



"L(B)^«A(B)- 



(F7) 



are the photon-number operators at the outputs of the 
polarizing beam-splitters PBS. The operator input- 
output relations of the polarization analyzers. 



aHA(B, cos^ACB) 



sint^A(B) 



sin0A(B), (F8) 

,COS0A(B), (F9) 



relate the operators 04^(3) with the field operators of the 
corresponding horizontal and vertical modes on the in- 
puts of the polarization analyzers, aH,„, and a^..^^, re- 

A(B) A(B} 

spectively. Finally, we utilize the input-output relation 



for the atmospheric channel, cf. Eq. (Al), and express 

the operators aHA(B) ^^'^ '^^a(b) field-operators 
a!" and dl" at the transmitter, 

"A(B) A(B) ' 



«Va(b) = 



^«HA(B,+V/1^^HA<3,- (FIO) 

, (Fll) 



where the fluctuating transmission coefficient T is equal 
for all four modes due to negligible depolarization in the 
Earth atmosphere and small time delay between entan- 
gled photons. 

A special part of the device is the 50:50 beam-splitter 
BS, cf. Fig. |Fl[ which randomly separates entangled pho- 
tons. The devise registers only events, when the photons 
left different output ports of the beam splitter. The other 
half of the events, when the photons come to the same po- 
larization analyzer, are not resolved by this device. For 
simplicity, we can consider an equivalent scenario, see 
Fig. 



F2 



Time-separated modes are replaced by spatial- 
ones. At each side the modes are split by two 50:50 beam 
splitters BS. The resolved events in the original device, 
cf . Fig. |F1[ correspond to the following events in Fig. |F2| 
(a) clicks at the reflection port in the side A and the 
transmission port in the side B, (b) clicks at the trans- 
mission port in the side A and the reflection port in the 
side B. The probability of both events is 2Pij^^i^ {6a, Ob), 
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which should be used instead of Pij^Ab (^a, ^b) with addi- 
tional 50% losses for the detection efficiency 77 caused by 
the bea m-s plitters BS. However, the factor 2 disappears 
in Eq. (F2). Hence, the main effect of the beam-splitter 



BS consists in additional 50% detection losses. 

The parametric down conversion source irradiates the 
states 



+00 



|PDC) = (coshx)"^ J2 V^i^tanh" x , (F12) 



n=0 

where x is the squeezing parameter, and 



(F13) 



1 " 

= E(-i) 



For n = 1, 



l*i> 



;J|(|1)hJ0)v. Io)hbIi)v, 

|0>Hjl)vJl)He |0>V3) 
1 



(F14) 



|H)aIV)b-|V>a|H)b 



is the BeU state, which for {e'^^\e'i\e'-l\9'^^) = 
f ' f ' ^) maximally violates the Bell inequality. It 
is convenient to solve the problem using the Glauber- 
Sudarshan P representation. The corresponding charac- 
teristic function for the state (F12) reads as 



*in (/3ha , , , /?vb) = exp 



X exp 



tanh^ X |/3v^ 1^ + tanh^ X |/3hb |^ - tanhx (/?va/3hb + P^Jn^) 
1 — tanh^ X 

tanh^ X |/3h^ I' + tanh^ X |/3vb |' + tanhx {Pn^Pv^ + PhJvJ 
1 — tanh^ X 



(F15) 



Similar to Appendix Bl we utilize Eqs. (F8|-(F11) and quantum-state input-output relations (B5). This results in the 



characteristic function at the outputs of the polarization analyzers. 



X 

^out (/3.^,/3«^,/3r3,/?«3) = J dTP(r)exp 



T2 tanh^ X{ I3tA + WnA + WtJ + Wn, 



1 — tanh X 



X exp 



X exp 



tanhx {Pt^Pt, + P*j;^ + Pnjn, + [3*, J* J sin (^b - ^a)" 
1 — tanh^ X 

tanhx {pTjn,, + /3*,/3;3 ~ PuJt^ ~ /3;,/3*J cos (0b - ^a) 
1 — tanh^ X 



In terms of the characteristic function, Eq. ( F5 1 is given by 

+00 

P....b(^a,^b) = / d«/3$out /3.b,/3h3) Kc{P^,)Kc{P^^)Ko{P,,)Ko{PJ^) , 



(F16) 



(F17) 



where d^(3 = <f P^^<f P^^d^ P^^<f P^^, 



ifo(/3) = — exp 

TTT] 



N 



(F18) 



KcW) ^ SiP) - KoW) . (F19) 
Substituting (FI6I in (F17I we obtain an explicit form for the probabilities of simultaneous clicks at both sides, 

exp {-2N) 



PiA,iB (0a,^b) = (l-tanh^x)' 



Co + 2Ci + a 



2cxp(-37V)\ /cxp(-4iV)\ 



Co + Ci / + \ Co I 



where (. . .) means averaging with respect to T, 

Co = {?7'r4 tanh^ x - [l + (r/T^ - l) tanh^ x] , 



(F20) 



(F21) 
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Ci ^ rjT^ (1 - rjT^) (l - tanh^ x) tanh^ x [v^T* tanh^ X - [l + (ryT^ - l) tanh^ x] ^} 



(F22) 



Ct^,Tb = C^Ra.Rb = V^T^ tanh^' x [l - tanh^ x] ' [(l - rjT^fta.nh' x - sin^ (Oa - ^b) 



(F23) 



Ta,Rb 



Ra,Te 



— jf'T^ tanh^ X [l ^ tanh^ x] 



77T^)^tanh^X-cos^' 



(F24) 



Equation (F20 ) can now be used in Eqs. (F2 |-(F4 ), and then inserted in expression for the BeU parameter (Fl I. Plots 
for the Be 



1 parameter as a function of the squeezing parameter x is given in Fig. 3. 
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FIG. F2. Equivalent scenario for the scheme presented in 
Fig. |F1[ The time separation is replaced by the spatial sepa- 
ration with two beam splitters BS at each side. The schemes 
(a) and (b) correspond to two kinds of events, which are 
detected by the original device. 



From the paper [7] we can extract some values of pa- 
rameters. Detection efficiency is 77 = 0.25 (6dB). Addi- 
tionally we have to include the beam splitter BS losses. 
The total efficiency of the receiver module is ry = 0.125 
(9 dB). The dark-count rate is 200 s~^ and the stray-light 
rate is 200 s"'^. Coincidences are integrated over a 1.25 ns 
time window. This gives the mean number of stray-light 
and dark counts of A^=5xl0~''. The mean losses of the 
atmospheric channel are (T^) =6.3x10"^ (32 dB). 



Appendix G: Quadrature and photon-number 
squeezing 

In this Appendix we consider some mathematical de- 
tails related to the post-selection procedure, which im- 
proves quadrature and photon-number squeezing. We 
will characterize the photon-number squeezing by the 



Mandel parameter. 



(: An2 :) 



(Gl) 



where h is the photon-number operator, : . . . : means 
normal ordering. Similarly, the quadrature squeezing is 
characterized by the normally-ordered variance 



of the quadrature 



X = a 



(G2) 



(G3) 



where a is the field annihilation operator. One can also 
use the value of squeezing in dB relatively to the vacuum 
noise. 



-loiogio (g + 1) 

for the photon-number squeezing, and 
-101ogio((: A^' 



(G4) 



(G5) 



for the quadrature squeezing. 

Input-output relation ( B8 1 for the normally-ordered 



moments can be applied for obtaining the input-output 
relations for squeezing factors (Gil and (G2), 



Qout ^ ji2^ Qii 



(Ml 



(: AX2 :)out = (r')(: AI^ :).^ + (AT2)(1) 



2 

in 7 



(G6) 



(G7) 



where (AT^) and (Ary^) are the variances of the trans- 
mission coefficient T and the transmission efficiency ry,, = 
T^, respectively, cf. also Refs. [TSl[2n]- The first terms in 
Eqs. (|G6|) and (G7) resemble the standard attenuation. 



The second terms are non-negative and increase with in- 
creasing and (nin)- They describe the effects of 
loss fluctuations. Hence the nonclassicality diminishes 
by two effects: the beam-positioning fluctuations and the 
brightness of the input light. 

An appropriate post-selection by the transmission co- 
efficient T can improve the nonclassicality of detected 
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fields. The corresponding procedure excludes the events 
with T < Tmin. The remaining events are described by 
the modified PDTC, 



J[0, else 



where 



(G8) 



In order to characterize the feasibility of the post- 
selection procedure, we can explicitly express the squeez- 



ing factors ([06]) and jGTj) [and post-selected PDTC ( [G8| ] 
via the exceedance J^. For this purpose we have to re- 
solve Eq. ( G9 ) with respect to T„iin and substitute it in 
post-selected PDTC ( C8 1 . For the particular case of d—0 



the corresponding relation, 



^ '^(-^min) 



(G9) 



is the value of exceedance ( E10[ ) at the point T—T^ 



Post-selected PDTC ( G8 1 can be used in Eqs 



( G7 1 for calculation of the improved values of t 



G6) and 



ne corre- 



sponding squeezing factors. 



TqCxp 



is obtained explicitly. 



2^2 
i?2 



2 ln(l--^) 



(GIO) 



